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Removing  Resonant  Oscillation  Signals 


from  Swallow  Float  Data 
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ABSTRACT 

Self-contained,  freely-drifting  Swallow  floats  capable  of  recording  very  low  frequency  (VLF)  ambient 
ocean  noise  are  under  development  at  the  Marine  Physical  Laboratory,  Scripps  Institution  of  Oceanogra¬ 
phy,  San  Diego,  California.  The  floats  are  ballasted  to  neutral  buoyancy  at  a  desired  depth,  whereupon 
they  record  the  components  of  particle  velocity,  from  which  sound  pressure  levels  may  be  derived.  A  high 
frequency  acoustic  mutual  interrogation  system  is  incorporated  to  provide  relative  float  positions,  permit¬ 
ting  the  combination  of  data  from  an  array  of  floats. 

During  an  experiment  conducted  between  16  and  18  September  1986  approximately  50  miles  west  of  San 
Diego,  twelve  Swallow  float  buoys  were  deployed  to  depths  of  1000  to  2000  meters  for  a  24  hour  period. 
Data  from  four  of  the  floats  were  found  to  contain  strong  signals  due  to  resonant  float  oscillation  excited  by 
their  internal  tape  recorder.  This  paper  develops  and  applies  a  method  of  processing  data  to  attenuate  the 
resonant  oscillation  signals.  The  method  can  be  applied  when  the  mean  square  resonant  oscillation  signal 
is  large  relative  to  the  background  noise  power. 

Consecutive  records  are  averaged  together  to  produce  an  estimate  of  the  float’s  resonant  oscillation  signal. 
The  resonant  oscillation  signal  estimate  is  then  correlated  with  each  record  to  yield  a  weighting  coefficient, 
and  the  weighted  signal  subtracted  from  that  record.  The  method  was  applied  to  three  of  the  four  1986 
deployment  floats,  in  which  it  significantly  attenuated  resonant  oscillation  signals  between  0.36  Hz  and  6.1 
Hz. 


Introduction 


Under  Office  of  Naval  Research  sponsorship,  the  Marine  Physical  Laboratory  has  designed,  fabri¬ 
cated  and  deployed  self-contained  Swallow  floats  which  record  very  low  frequency  (VLF)  ambient  ocean 
noise  over  extended  periods  of  time.  In  operation,  the  floats  are  ballasted  to  neutral  buoyancy  at  a  desired 
depth  or  deployed  to  the  ocean  bottom.  They  operate  autonomously  to  measure  and  record  the  components 
of  particle  velocity  in  the  1  -  20  Hz  band,  from  which  sound  pressure  levels  may  be  determined.  The  floats 
periodically  generate  and  receive  8  kHz  acoustic  pulses  which  are  used  to  determine  their  relative  posi¬ 
tions.  The  deployment  of  several  floats  thus  forms  a  freely  drifting  array  of  sensors. 

The  Swallow  float  design  minimizes  self-noise  which  can  limit  accurate  ambient  ocean  noise  meas¬ 
urements.  Floating  freely,  they  are  not  subject  to  flow  noise  or  tether  strumming.  They  measure  velocity 
and  are  thus  insensitive  to  variations  in  local  pressure.  Ambient  ocean  noise  levels  derived  from  single 
float  measurements  have  been  shown  to  agree  closely  with  omni-directional  hydrophone  measurements 
made  under  similar  conditions.1  Operation  of  the  floats  has  been  described  in  earlier  reports. '•2- 3  Starting 
with  the  1986  deployment  data,  the  velocity  signal  sampling  rate  was  been  increased  to  SO  Hz  and  the 
anti-aliasing  filter  cutoff-frequency  increased  to  about  21  Hz.  Record  duration  was  decreased  from  90  to 
45  seconds. 

MPL  has  deployed  Swallow  floats  annually  since  1982.  Between  16  and  19  September  1986,  12 
floats  were  deployed  approximately  SO  miles  west  of  San  Diego,  California  at  32°3S'  N,  118°10'  W. 
Three  floats  were  located  on  the  ocean  bottom  to  serve  as  references  for  the  acoustic  positioning  system. 
Five  floats  developed  problems  and  yielded  little  data.  The  remaining  4  deployed  to  about  1000  meters 
depth  where  they  acquired  approximately  IS  hours  of  acoustic  data.  MPL  Technical  Memorandum  391 
contains  a  preliminary  analysis  of  data  from  this  deployment 

Velocity  data  from  the  floats  deployed  to  1000  meters  depth  in  1986  have  been  found  to  contain 
strong  sinusoidal  signals  which  are  excited  at  the  start  of  each  record  and  decay  exponentially.  These  sig¬ 
nals  are  attributed  to  resonant  float  oscillation  apparently  excited  by  the  internal  tape  recorder  when  it 
operates  between  records.  In  some  records,  the  resonant  oscillation  causes  velocity  data  to  be  clipped  for 
up  to  8  seconds  at  the  beginning,  obscuring  the  ambient  ocean  noise  signal  and  reducing  the  dynamic  range 
with  which  the  ambient  ocean  noise  signal  is  measured  during  the  remainder  of  the  record. 

Effort  is  currently  underway  at  MPL  to  modify  the  Swallow  float  hardware  to  eliminate  resonant 
oscillation  in  future  deployments.  It  is  also  desirable  to  remove  these  signals  from  the  1986  data  set  to  per¬ 
mit  its  further  use.  This  report  derives  a  method  of  removing  or  significantly  reducing  resonant  oscillation, 
which  may  be  applied  if  certain  conditions  are  met.  The  method  is  shown  to  be  effective  on  the  1986 
deployment  data  set. 

A  separate  investigation  focusing  on  occurrences  of  resonant  oscillation  at  frequencies  below  0.6  Hz 
is  ongoing.  Resonant  oscillation  at  these  frequencies  appears  to  consist  of  an  irregular  rocking  or  wobble 
which  is  strongly  excited  as  the  floats  descend,  weakly  but  consistently  excited  by  the  tape  recorder  at  the 
start  of  each  record,  and  occasionally  excited  by  strong  acoustic  pulses  which  it  receives.  MPL  Technical 
Memorandum  394,  Low  Frequency  Oscillations  in  1986  Swallow  Float  Data ,  will  probably  be  issued  in 
June  of  this  year. 
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1.  Swallow  Float  System  Description 

The  Swallow  float  is  a  17  inch  diameter  glass  sphere  containing  three  geophones  as  directional  velo¬ 
city  transducers,  a  compass  for  detecting  float  heading,  an  acoustic  transponder  for  positioning,  a  solid  state 
memory  data  buffer,  a  digital  cartridge  tape  recorder,  and  an  acoustically  actuated  ballast  release.  Signal 
conditioning  is  depicted  in  Figure  1.1  and  hardware  configuration  is  shown  is  Figure  1.2. 

The  three  orthogonally-oriented  geophones  are  rigidly  mounted  to  the  inside  top  of  the  glass  sphere. 
Each  transducer  generates  a  signal  proportional  to  the  component  of  the  float’s  velocity  which  is  along  the 
transducer’s  axis. 

The  geophone  signal  undergoes  95  dB  of  amplification  before  entering  an  automatic  gain  control 
(AGC)  circuit  which  provides  up  to  26  dB  of  gain  to  make  full  use  of  the  dynamic  range  of  (he  analog-to 
digital  converter  (ADC).  AGC  level  is  controlled  by  the  float's  microprocessor  depending  upon  how  much 
clipping  occurs  at  the  ADC. 

Before  sampling  at  50  Hz,  the  signal  is  low-pass  filtered  to  prevent  aliasing.  The  filter  has  a  21  Hz 
comer  frequency  and  is  down  20  dB  at  25  Hz.  The  ADC  is  an  8-bit  device  yielding  48  dB  of  dynamic 
range.  Clipping  occurs  when  the  input  rises  above  +  2.5  or  below  -  2.5  volts. 

The  digitized  output  is  buffered  for  44  seconds  before  being  stored  on  tape  because  the  tape  drive 
requires  approximately  75  watts  of  power  and  cannot  be  powered  on  continuously.  During  a  one  second 
interval,  the  tape  drive  is  powered  up,  the  buffer  contents  are  written  to  tape,  and  the  tape  drive  is  powered 
down.  Tape  recorder  capacity  is  limited  to  about  2000  45-second  records,  so  that  submergence  time  is  lim¬ 
ited  to  about  25  hours. 


Figure  1.1:  Geophone  Signal  Conditioning 


2.  The  Resonant  Oscillation  Problem 


Twelve  Swallow  floats  were  deployed  during  September  1986  about  50  miles  west  of  San  Diego,  Ca. 
in  approximately  2000  meters  of  water.  Three  were  located  on  the  ocean  bottom  to  serve  as  references  for 
the  acoustic  positioning  system  and  the  remaining  nine  were  deployed  to  about  1000  meters  depth  to  col¬ 
lect  velocity  data  Four  of  the  nine  surfaced  prematurely  and  collected  little  data.  A  fifth  float  developed  a 
ground  in  one  geophone  channel.  The  four  remaining  floats  deployed  successfully  to  approximately  1000 
meters  and  collected  about  15  hours  of  useful  velocity  data.  A  preliminary  analysis  of  these  data  is  con¬ 
tained  in  MPL  Technical  Memorandum  391,  September  1986  Trip  Report* 

The  1986  deployment  velocity  data  were  found  to  contain  sinusoidal  signals  which  are  strong  at  the 
start  of  each  record  and  decay  exponentially  thereafter.  Some  examples  of  geophone  output  with  95  dB 
gain  are  shown  in  Figures  2.1  through  2.8.  Each  Figure  contains  one  axis  of  data  from  twelve  45-second 
records.  Acoustic  positioning  pulses  may  be  seen  10  seconds  into  records  1203,  1611,  1205,  1613,  1209. 
1605, 1210,  and  1606  in  Figures  2.1  through  2.8,  respectively. 


Table  2.1:  Key  to  Figures  2.1  through  2.8 

Figure 

Float 

Records 

2.1 

3 

1200-1211 

2.2 

3 

1600-1611 

2.3 

5 

1200-1211 

2.4 

5 

1600-1611 

2.5 

9 

1200-1211 

2.6 

9 

1600-1611 

— 

10 

1200-1211 

2.8 

10 

1600-1611 

Floats  3,  9,  and  10  contain  a  =4  Hz  signal  which  is  clipped  for  3  to  8  seconds  at  the  start  of  each 
record  and  decays  away  after  that  Another  lower  amplitude  sinusoid  with  a  period  of  about  2.6  seconds  is 
visible  in  many  records  in  all  4  floats.  The  decaying  sinusoids  appear  quite  similar  from  record  to  record. 
Their  initial  phase  does  not  change  and  their  amplitude  changes  very  little.  For  these  reasons,  they  are 
thought  to  represent  float  resonances  excited  by  one  or  more  impulses  which  are  nearly  identical  from 
record  to  record.  The  impulses  are  thought  to  be  mechanical  rather  than  acoustic.  This  is  because  the 
floats  do  not  generate  acoustic  energy  near  the  start  of  each  record,  and  because  mechanical  impulses  are 
quite  possibly  generated  by  the  float’s  internal  tape  recorder  when  it  is  powered  up  and  down  at  the  end  of 
each  record. 

Another  indication  that  the  tape  recorder  is  responsible  for  the  oscillations  is  that  the  initial  amplitude 
and  phase  of  the  oscillation  depends  on  the  direction  in  which  the  tape  recorder  is  writing.  The  tape 
recorder  writes  in  a  serpentine  manner,  i.e.  it  traverses  the  entire  tape  in  one  direction  writing  the  first  track, 
reverses  direction  and  writes  the  second  track,  again  reverses  direction  and  writes  the  third  track,  and  so  on. 
Each  track  is  about  500  records  long,  with  direction  reversals  occurring  near  records  500,  1000  and  1500. 
An  example  of  the  effect  which  tape  recorder  direction  has  on  the  resonant  oscillation  can  be  seen  in  data 
from  float  3.  X-axis  clipping  near  record  1200  (track  3)  in  Figure  2.1  has  an  initial  phase  which  is  about  Jt 
radians  different  from  (hat  near  record  1600  (track  4)  in  Figure  2.2.  Similar  indications  can  be  found  in  the 
other  floats. 

The  resonant  float  oscillation  causes  at  least  two  problems.  First,  it  can  obscure  the  ambient  noise 
signal.  Certainly  the  ambient  noise  signal  is  obscured  at  the  start  of  (he  record  in  floats  3,  9,  and  10.  How 
much  of  the  record  is  affected  by  the  oscillations  it  is  not  readily  evident.  The  0.4  Hz  oscillation  appears  to 
extend  to  the  end  of  some  records. 
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Second,  resonant  oscillation  can  dominate  the  AGC  level  and  reduce  the  dynamic  response  applied 
to  the  ambient  noise  signal.  To  understand  why  this  can  occur  requires  some  understanding  of  how  the 
AGC  level  is  set.  The  AGC  level  is  controlled  by  a  microprocessor  which  keeps  track  of  the  number  of 
clipped  samples  during  each  record.  The  AGC  is  stepped  down  between  records  when  there  is  too  much 
clipping  and  is  stepped  up  when  there  is  too  little.  The  resonant  oscillation  can  raise  the  AGC  level  by 
causing  clipping.  When  the  AGC  level  is  set  too  high,  too  few  bits  are  available  to  quantize  the  ambient 
ocean  noise  signal.  This  reduces  the  dynamic  range  with  which  the  noise  is  measured.  The  microprocessor 
can  be  programmed  to  wait  after  the  start  of  a  record  before  looking  for  clipped  samples.  This  works  well 
enough  but  requires  knowing  how  long  the  oscillation-induced  clipping  will  last. 

Fundamentally  there  are  two  ways  of  modifying  the  Swallow  float  design  to  eliminate  resonant  oscil¬ 
lation  in  the  velocity  data:  remove  the  excitation  or  change  the  float’s  response  to  the  excitation.  Eliminat¬ 
ing  the  excitation  is  impractical  because  it  would  be  expensive  to  redesign  the  tape  recorder  or  would 
require  long,  very  soft  "springs"  to  decouple  it  from  the  rest  of  the  float  down  to  1  Hz.  Changing  the 
response  is  easier  done  and  effort  is  ongoing  at  MPL  to  dampen  the  resonant  response,  move  it  out  of  band, 
or  both.  One  design  was  tested  during  a  February  1987  sea  trip  but  proved  unsatisfactory.  Several  other 
designs  will  be  tested  during  an  April  1987  sea  trip. 

There  is  also  sufficient  value  in  the  September  1986  data  set  to  warrant  exploring  a  means  of  pro¬ 
cessing  that  data  to  remove  the  resonant  oscillation.  The  intent  with  the  1986  data  is  to  use  the  acoustic 
positioning  system  pulses  to  determine  the  floats  positions  while  they  were  deployed  and  beamform  using 
data  from  the  four  floats  deployed  to  1000  meters.  The  remainder  of  this  report  deals  with  removing  the 
resonant  oscillation  from  the  1986  data. 


3.  Theoretical  Solution 


The  j-  Swallow  float  data  record  x,  (n )  is  composed  of  an  ambient  noise  signal  sj  (n )  and  a  resonant 
oscillation  signal  ry(n).  The  resonant  oscillation  r;(n)  consists  of  a  fixed  but  unknown  sequence  r(n), 
multiplied  by  a  gain  factor  ay  which  may  vary  from  record  to  record.  See  Figure  3.1.  The  data  record 
is  therefore 

Xy(rt)  =  ayr(/j)  +  jy(fl).  (3.1) 


We  define  r(n)  such  that  the  average  value  of  fly  is  1.  It  can  be  seen  from  Figures  2.1  through  2.8  that  the 
value  of  fly  changes  little  from  record  to  record.  Therefore  we  assume  that  for  M  greater  than  10, 

a  =  aJ  =  l  (M  >  10)  (3.2) 

M  y= 1 


Figure  3.1:  The  geophone  signal  x(n)  is  composed  of 
an  ambient  noise  signal  s(n)  and  a  resonant  oscillation  signal  r(n). 

A  method  is  sought  by  which  the  resonant  oscillation  ayr(n)  may  be  removed  from  each  record 
Xj(n)  leaving  behind  the  noise,  s,(n ),  which  is  the  desired  signal.  The  approach  is  to  estimate  r(n)  and  ay 
and  substitute  them  into  a  rearranged  equation  (3.1): 

*/(n)  =  Xj(n)  -  dj  f(n)  (3.3) 


Estimating  r(n) 

Under  certain  conditions,  a  good  estimate  of  r  (n )  may  be  obtained  by  coherently  averaging  data 
from  a  large  number  of  records.  The  M  record  average  of  x,  (n )  is 
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1  *  1  *  r  1 

c(")s77Z;ty(n)=  77  Z  a>T(*  )  +  *>(«) 
My,,  M  jm !  L  J 


=r(n)  TrZ  a/  +  7 jZ  Sj(n)ma-r(n)  +  s(n) 

M  J.l  M  jml 


Equation  (3.4)  tells  us  that  if  M  is  large  enough  for  a  to  be  approximately  equal  to  1 ,  and  if  I  r  (n )  I  is  much 
larger  than  I s(n )  I ,  then  r(n)  will  be  approximately  to  x(n ). 

If  lr(n)l  »  I  J(n)l,  (3.5) 


1  M 

then  r(n)=  —  £  Xj(n)  =  r(n) 

M  *-i 


(M  >  10) 


We  can  ensure  that  equation  (3.5)  is  satisfied  in  the  mean  square  sense.  Squaring  both  sides  and  averaging 
over  all  values  of  n  yields 


T  I  r\n)  »  yZ  s\n) 


Using  equations  (3.5a)  and  (3.6),  we  can  say  that  if  the  mean  square  value  of  r(n)  is  much  greater  than  the 
mean  square  value  of  s(n ),  then  the  M  record  average  provides  a  good  estimate  of  r  (n ).  We  need  to  know 
the  mean  square  values  of  r  (n )  and  T(n ). 

The  sequence  r  (n )  consists  of  fixed  but  unknown  values  which  may  be  positive,  negative  or  zero. 
No  assumption  is  made  about  its  mean  value 

1  L 

Mr  =  hm  —  £  r(n)  (3.7a) 

L-¥  ••  L  * 


except  that  it  is  well  approximated  by  the  estimate  \ir 


1  L 

Mr  =  7  E  '■('*)  =  Mr 
L  »=i 


(L  =  2250) 


when  L  =  2250,  the  number  of  points  in  a  record.  There  is  no  a  priori  assumption  about  the  mean  square 
value  ofr(n). 


Yr2  =  lim  7  E  r2(n) 

L~*  “  L  »=i 


except  that  it  is  well  approximated  by  the  estimate 


P 


|L- 

* 


>  a 
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#=T  I'V)  =  Yr  (L  =  2250) 

L  »=i 


The  sequence  Sj(n)  is  one  realization  of  a  random  process  which  we  assume  to  be  stationary  over  M 
record  intervals.  We  assume  that  st  (n )  is  mean  zero 

1  M 

M.  =  —  2  *,(*)  =  0  (3.9a) 


and  has  mean  square  value 


t  M 

1/2-  1  i  m  -  V  e  2/ 


The  stationarity  assumption  also  allows  us  to  equate  the  time  and  ensemble  averages  of  s,  (n ),  and  we 
assume  that  they  are  well  approximated  by  the  L  =  2250  estimate. 


=  7-  E  */(«)  =  0 

^  «=i 


(L  =  2250) 


=  (L  =  2250) 


We  also  assume  that  the  noise  is  uncorrelated  from  record  lo  record 


1  ^  f  0* 

^*(0)=  lun  Sj(n)sk(n)  =  j  2  .  . 


and  that  the  estimate  of  the  cross  correlation  (with  L  =  2250)  is 


(3.10a) 


(3.10b) 


(3.11a) 


-  i  l  r  o  j  *k 

f?x,x.(0)=  -£  r;(n)  s*(n)  =  •  ^  ^ 


Lastly,  we  assume  that  r  (n )  and  s,  (n )  are  uncorrelated  for  all  j 


(L  =  2250) 


(3.11b) 


fln.(0)  =  ,lim  —  £  r(n)sk(n)  =  0 


(for  all  j) 


(3.12a) 


i 


m'm 
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/?„.( 0)=  fZ  r(n)sk(n)=0 


(L  =  2250;  for  all  j) 


(3.12b) 


Now,  the  estimate  of  the  mean  square  value  of  s(n )  is 

y?=j-Zs2(») 

u  «= i 


M 


M 


;=i 


1  M  M 

772  I  2  *,(»)**(«) 

M  j~\ 


,  Af  W  1 

772 II I 

M  /= 1  *=l  i- 


L 


I 


j/(n)sJk(/t) 


=  -L  v72=  A 

A/2  M 


(3.13) 


where  equation  (3.11b)  has  been  used.  Equation  (3.13)  says  that  the  mean  square  value  of  s(n)  is  1/M 
limes  the  mean  square  value  of  s  ( n ).  An  estimate  for  the  mean  square  value  of  s  (n ),  yj,  may  obtained 
from  the  mean  and  mean  square  values  of  xt  (n )  as  follows.  The  mean  value  of  the  j-  data  record  is 


fV 


L  ^ 

u  *=1 


E  Xj(n)  =  a j  —  £  '■(«)+  T  Z  *>(*)  =  <»>  1L- 


(3.14) 


where  equations  (3.7b)  and  (3. 10a)  have  been  used.  The  mean  square  value  of  x,  (n )  is 

Yi,  =  f  I  *A*> 


I  L 

—  £  a,r{n)  +  Sj(n) 
L  *=1  L 


2  1 

=  a,  - 


I  '2(")+f  I*A«) 


*=l 


2  * 2  A2 
=  «/  Yr  +  YT 


(3.15) 


using  equation  (3.12b).  Now  solve  equation  (3.14)  for  ar  square  the  result  and  substitute  it  into  equation 
(3.15)  to  get 


MPL-U-21/87 


*2  A2  Yr  A? 

— +7? 

Mr 


There  are  M  such  equations  and  they  can  be  written  in  matrix  form  as  follows 


vi  i  f— i 


If  we  identify 


so  that  equation  (3.17)  may  be  written  as  Z  =  HX ,  then  the  least  squares  solution  for  X  is 


=  \hth  x  H7 


**  /-i 

v2  _  « 

w  least  square  £  p.  ‘  A/ 

/=1 


M  M  I  _1  r  M  1 
Z^Ya’ 

/•I  y=i 


^  .721  I"  £  2  -2I 


M  ~  Z  4a,  Z  4a,  Ya, 
/=!  /  =  I 


M  M 

m  y  u?  -  y.  u.2 


M  M 

V'  *2  a-  *4 


~  Z  4a,  L4a4,  Z  Ya, 
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The  least  squares  estimate  for  the  mean  square  value  of  the  noise  is 


MM  MM 

V  V'  J2  V  ..2  v  ,.2  ..2 


2  K  I  -  Z  K  2  K  % 

:i  /“i  /-I 

&  = - u - n: — n — 


M  M 

>•>  .  /-1 


Substituting  (3.19)  into  (3.13)  provides  the  estimate  for  the  mean  square  value  of  the  J(n).  An  esti¬ 
mate  for  Yr  may  be  obtained  from  the  mean  square  value  of  the  M  record  average,  x(n).  Using  (3.4)  and 
(3.12b),  we  have 


i?>  r  i  [*<»>] 


-rz[* 

L  nml  L 


r(n)  +  s(n) 


1  L 

—2  *  T' 


=  *T  l'2<*>+f 


=  a2y?  +  % 


If  M  is  larger  than  10,  we  can  set  a  -  1  in  equation  (3.20).  Solving  (3.20)  for  yr2  yields 


Y?s  H  ~ir  M>1° 


Recognizing  that  the  condition  (3.5a)  is 


Y?  »  7?  • 


we  can  use  equations  (3.13),  (3.19)  and  (3.21)  to  determine  if  (3.5b)  is  satisfied  for  a  particular  value  of  M. 
If  (3.5b)  is  satisfied,  then  we  are  justified  in  using  that  value  of  M  in  equation  (3.6)  to  obtain  r  (n ). 


R 


EX 
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Estimating  a 

An  estimate  for  the  value  of  a  in  the  jth  record  may  be  obtained  by  correlating  xj  (n )  with  t  (n ): 


Rx,r( 0)  =  Y  E  xj(n)r(n) 

L  «=i 


=  7!  [ayrOO  +  jjOt)]  r(n) 
L  1  J 


:a/  T  £  r(n)r(/i)+  y  £  r(n)sj(n) 

L  11=1  L  <i=i 


*2 

=  a,Y, 


where  the  assumption  (3.12)  that  r  (n )  and  Sj(n )  are  uncorrelated  has  been  used.  Equation  (3.22)  is  valid  if 
r(n)  is  a  good  estimate  of  r(n),  which  is  to  say  that  the  condition  (3.5b)  is  satisfied.  Solving  (3.22)  for  a, 
yields 


Rv( 0) 


Equation  (3.23)  is  evaluated  by  calculating  y?  using  (3.21)  and  RXir(Q)  using  (3.22).  Substituting 
equations  (3.6)  and  (3.23)  into  (3.3)  yields  the  estimate  of  the  ambient  ocean  noise  in  the  jth  record 


/MO) 


s(n)  =  Xj(n)-djf(n)  =  Xj(n ) - —  x(n) 

Y r 


W*. 
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4.  Application 

In  the  previous  Section,  a  method  has  been  developed  for  removing  resonant  oscillation  signal  from 
Swallow  float  data  records.  In  this  Section,  the  method  is  applied  to  September  1986  deployment  data. 


Estimating  r(n) 

A  necessary  first  step  is  to  clarify  the  definition  of  r(n).  In  Section  3  it  was  postulated  that  the 
resonant  oscillation  r(n)  be  a  fixed  but  unknown  sequence  which  is  different  for  each  float.  In  Section  2, 
however,  it  was  pointed  out  that  the  initial  phase  of  the  resonant  oscillation  depends  upon  the  direction  in 
which  the  tape  recorder  is  writing.  In  fact,  r(n)  is  thought  to  change,  but  change  very  little  over  reason¬ 
ably  long  intervals  which  do  not  include  changes  in  tape  recorder  direction,  r  (n )  is  characterized  as  fol¬ 
lows: 

•  the  resonant  oscillation  r(n)  results  from  impulses  generated  by  the  tape  recorder  when  it  starts 
and  stops  between  records; 

•  the  magnitude  and  direction  of  the  impulses  generated  by  the  tape  recorder  depends  upon  the  direc¬ 
tion  in  which  the  tape  recorder  is  writing,  and  also  depends  upon  the  relative  amount  of  tape  which  is 
on  each  reel  of  the  recorder; 

•  the  magnitude  and  initial  phase  of  the  resonant  oscillation  depend  upon  tape  travel  direction  and 
may  vary  during  the  interval  between  tape  recorder  direction  changes. 

Thus  r(n  )  is  expected  to  change  considerably  when  the  tape  recorder  changes  direction,  and  to 
change  slowly  as  the  tape  moves  from  one  reel  to  the  other.  The  implication  of  this  redefinition  is  that  M 
record  averages  which  are  meant  to  approximate  r  (a)  cannot  span  tape  recorder  direction  changes  or  be  so 
long  that  r  (n )  changes  appreciably  over  the  average. 

Tape  recorder  direction  changes  have  been  shown  to  occur  at  records  968  and  1456  in  float  9  and  at 
records  977  and  1467  in  float  10.5  We  can  assume  that  tape  recorder  direction  changes  occur  in  floats  3  and 
5  at  roughly  the  same  points  and  avoid  averaging  over  those  records.  Once  we  have  determined  the  value 
of  M  which  satisfies  equation  (3.5b)  for  each  float,  we  can  ensure  that  r(n)  does  not  change  appreciably 
over  the  M  record  average  by  comparing  adjacent  M  record  averages. 

We  now  proceed  to  process  1986  deployment  data  from  floats  3,  5, 9,  and  10  and  start  by  finding  the 
values  of  M  which  result  in  the  mean  square  value  of  the  M  record  averages  being  much  greater  than  1/M 
times  the  mean  square  values  of  the  noise.  This  is  the  condition  expressed  in  equation  (3.5b). 

One  must  establish  how  large  the  left  hand  side  of  (3.5b)  is  to  be  in  order  that  the  condition  be 
deemed  satisfied.  We  started  with  the  condition  (3.5)  which  required  only  that  I  r  (n )  I  be  much  larger  than 
lj(«)l.  It  is  now  stipulated  that  I  r  (n )  I  must  be  at  least  a  factor  of  5  larger  than  ls(n)l. 


ls(n)l 


(4.1) 


Since  equation  (3.5b)  is  the  mean  square  of  (3.5),  equation  (4.1)  requires  that  y2  be  at  least  a  factor  of  25 
larger  than  y|. 


The  problem  is  now  to  find  the  value  of  M  for  each  float  which  satisfies  equation  (4.2).  Tables  4.1 
through  4.4  contain  the  results  of  evaluating  equation  (4.2)  for  various  values  of  M  for  floats  3,  5,  9,  and 
10.  The  right-most  column  is  each  Table  is  the  value  of  the  left  hand  side  of  equation  (4.2). 

According  to  Table  4.1,  setting  M  equal  to  60  satisfies  (4.2)  for  all  three  axes  of  float  3.  Table  4.2 
shows  that  using  M  =  400  satisfies  (4.2)  in  the  x  direction  only  for  records  1001-1400  of  float  5.  We  can¬ 
not  use  a  larger  value  of  M  without  running  the  risk  of  including  records  which  contain  tape  recorder  direc¬ 
tion  changes.  Thus  we  cannot  average  enough  float  5  records  to  satisfy  equation  (4.2).  Tables  4.3  and  4.4 
show  that  a  value  of  M  =  40  satisfies  (4.2)  for  all  axes  from  floats  9  and  10. 

From  Tables  4.1,  4.3  and  4.4  we  know  how  many  records  must  be  averaged  together  to  get  a  good 
estimate  of  r(n)  for  floats  3,  9,  and  10.  However,  before  going  on  to  process  float  3,  9,  and  10  data,  we 
must  also  ensure  that  r(n)  does  not  change  significantly  over  the  40  or  60  record  intervals  which  do  not 
include  tape  recorder  direction  changes.  To  do  this  we  determine  the  similarity  of  consecutive  M  record 
averages  as  indicated  by  their  covariance.  The  covariance  of  the  two  sequences  normalized  by  the  product 
of  their  variances  produces  their  correlation  coefficient  p.6  The  correlation  coefficient  of  two  M  record 
averages,  r*  (n)  and  rk(n),  is 


<V.(Q) 

Kjfj(oycM(oj 


where: 


O,*(0)  =  Y  I  r'M)~  14,]  ['*(«)  - 14.] 

<V,(0>  =  l  I  [r> >  - 14,]  [rj (* )  - 14,]  . 

Qr.(0)  =  £  Z  ['*(«)  *  14.]  ['*(n )  - 14, ] 


1  L 

14,  =  Itm  —  £  r"(/i)  and 

L  -4  ••  L  _  _  i 


1  L 

It,.  =  lim  —  Z  '>(")• 
L~*  ~  L  «=! 


When  L  =  2250  is  used  in  equations  (4.1)  through  (4.6),  the  estimates  p,*,  Cr  r-(0),  Cf  f  { 0),  Cr._^(0),  pr/, 
and  (ip,  are  obtained. 
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The  correlation  coefficient  p/t  is  normalized  to  take  on  values  between  -1.0  and  1.0.  A  value  of  p,* 
close  to  1.0  indicates  that  the  two  sequences  are  very  similar.  A  value  close  to  0  indicates  that  they  are 
unrelated,  and  a  value  close  to  -1.0  indicates  that  they  are  inversely  related. 


Figure  4.1  shows  the  correlation  coefficients  for  a  series  of  adjacent  60  record  averages  from  float  3. 
The  correlation  coefficients  are  very  close  to  1  everywhere  except  in  the  vicinity  of  the  tape  recorder  direc¬ 
tion  change,  which  apparently  occurs  somewhere  between  record  1360  and  1420.  Figure  4.1  indicates  that 
r  ( n )  changes  little  over  60  record  intervals  which  do  not  include  a  change  in  tape  recorder  direction. 


Figures  4.2  and  4.3  show  the  correlation  coefficients  for  series  of  adjacent  40  record  averages  from 
floats  9  and  10,  respectively.  Again  the  correlation  coefficients  remain  very  close  to  1  everywhere  except 
in  the  vicinity  of  the  tape  recorder  changes  which  apparently  occur  between  records  1400  and  1440  in  float 
9  and  records  1440  and  1480  in  float  10.  Figures  4.2  and  4.3  indicate  that  r  (n )  in  floats  9  and  10  changes 
little  over  40  record  intervals  which  do  not  include  a  change  in  tape  recorder  direction. 


We  now  have  full  confidence  in  using  a  60  record  average  as  the  estimate  of  r(n)  in  float  3  and  40 
record  averages  as  estimates  of  r  (n )  in  floats  9  and  10.  Before  moving  on  to  estimating  a,  and  processing 
the  data,  it  is  interesting  to  compare  estimates  of  the  mean  square  noise,  or  noise  power,  from  the  four 
floats  during  the  same  intervals,  summarized  in  Table  4.S. 


There  is  general  agreement  among  measurements  made  by  the  floats,  although  there  are  some  differ¬ 
ences.  The  floats  were  separated  by  1900  to  3000  meters  during  the  two  intervals  and  their  local  sound 
fields  may  not  be  the  same.  Also,  there  is  probably  variation  among  the  sensor  frequency  responses.  X  and 
y  direction  power  estimates  should  generally  agree,  as  they  do  in  floats  3,  5,  and  9  during  records  1201- 
1240  and  float  5  during  records  1601-1640,  since  all  four  floats  were  not  oriented  in  the  same  direction 
when  these  measurements  were  made.  When  they  do  not  agree,  as  in  float  10  during  records  1201-1240 
and  floats  3,  9,  and  10  during  record  1601-1640,  the  y  direction  levels  are  higher,  perhaps  due  to  the 
influence  of  the  resonant  oscillation  signal. 
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have  comparable  power  in  their  resonant  oscillation. 


Table  4.6  •  Resonant  Oscillation  Power  Estimates 


Float 

X 

Axis 

y 

z 

records  1201-1240 

float  3 

0.0223 

0.0606 

0.0044 

float  5 

0.0036 

0.0022 

0.0004 

float  9 

0.0255 

0.0444 

0.0084 

float  10 

0.0567 

0.1775 

0.0251 

records  1601-1640 

float  3 

0.0260 

0.0622 

0.0033 

floats 

0.0036 

0.0027 

0.0004 

float  9 

0.0273 

0.0948 

0.0100 

float  10 

0.0563 

0.1985 

0.0291 

Estimating  a 

Now  we  can  calculate  values  for  d,  using  equation  (3.23).  Values  of  a,  for  records  1201-1240  and 
1601-1640  from  floats  3, 9,  and  10  are  plotted  in  Figures  4.4  through  4.6.  The  mean  and  standard  deviation 
of  the  a;  values  are  given  in  Table  4.7.  The  mean  value  of  a;  is  very  close  to  1 .0  for  all  three  floats  and  the 
standard  deviation  never  exceeds  0.1.  Thus  is  appears  we  are  justified  in  our  assumption  (3.2). 


Table  4.7  •  Statistics  of  the  scaler  multiplier  a 


float 

axis 

records  1201-1240 

records  1601-1640 

3 

X 

Ml 

— 

0.090 

y 

Z 

0.9970 

9 

X 

1.0000 

0.037 

1.0000 

0.029 

y 

1.0000 

0.047 

1.0000 

0.028 

Z 

1.0000 

0.063 

1.0000 

0.034 

10 

X 

1.0000 

0.066 

1.0000 

0.057 

y 

1.0000 

0.025 

1.0000 

0.023 

Z 

1.0000 

0.051 

1.0000 

0.020 
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Processing  Results 

Data  from  floats  3, 9,  and  10  have  been  processed  using  M=40  or  60  and  equation  (3.24).  Results  are 
contained  in  Figures  4.7  through  4.12. 

Figures  4.7a  through  4.7h  pertain  to  float  3.  Figure  4.7a  shows  the  x,  y,  and  z  axis  geophone  time 
series  from  record  1207  before  any  processing.  These  same  time  series  were  shown  as  the  8th  traces  from 
the  top  in  Figures  2.1a,  2.1b,  and  2.1c.  The  time  series  consist  of  signal  due  to  resonant  float  oscillation, 
which  is  clearly  visible  in  the  first  3  to  10  seconds  of  the  record,  and  signal  due  to  ambient  ocean  noise. 

Figure  4.7b  contains  the  power  spectrum  of  the  time  series  shown  in  Figure  4.7a  The  first  5  seconds 
of  the  record  were  not  used  in  the  spectral  estimate  because  data  are  clipped.  A  512  point  (20.48  second) 
Fourier  transform  was  used,  and  consecutive  segments  which  were  50%  overlapped.  Segments  were  win¬ 
dowed  using  a  Kaiser-Bessel  window  with  a  =  2.5.  The  seven  spectral  estimates  thus  obtained  from  each 
axis  of  the  record  were  incoherently  averaged  to  produce  one  spectrum  per  axis.  The  acoustic  pressure 
spectral  estimate  shown  in  the  top  panel  was  obtained  by  power  summing  spectra  from  each  axis.  It  has 
been  shown  p  viously  that  acoustic  pressure  spectra  may  be  derived  from  the  power  sum  of  directional 
velocity  spectra. 1  Peaks  in  the  geophone  and  acoustic  pressure  power  spectra  may  correspond  to  signals 
present  in  either  the  resonant  oscillation  or  the  ambient  ocean  noise. 

Figure  4.7c  shows  the  lime  series  from  a  60  record  average  beginning  with  record  1201.  This  is  our 
estimate  of  the  resonant  oscillation,  f{n ).  The  first  3  to  10  seconds  of  r(n )  closely  resemble  record  1207, 
indicating  that  the  first  part  of  record  1207  is  strongly  affected  by  the  resonant  oscillation.  Figure  4.7c 
shows  that  the  resonant  oscillation  signal  continues  to  be  present  throughout  the  record,  although  its  level  is 
low  in  the  latter  part 

Figure  4.7d  shows  the  power  spectra  of  the  time  series  shown  in  Figure  4.7c.  The  levels  of  spectral 
peaks  are  summarized  in  Table  4.8.  Only  the  peaks  at  2  Hz  and  multiples  of  2  Hz  appear  harmonically 
related.  The  8  Hz  peaks  in  the  x  and  y  direction  are  2  Hz  harmonics  but  they  are  dominated  by  another  sig¬ 
nal,  as  evidenced  by  their  higher  levels. 

Figure  4.7e  compares  the  power  spectra  of  record  1207  with  that  of  r(n).  The  background  level  of 
r'(n  )’s  spectra  is  about  18  dB  below  that  of  record  1207.  This  difference  makes  sense  in  light  of  equation 
(3.13),  which  says  that  the  power  in  s(n)  is  1/M  times  that  in  s(n),  and 

10  log  ^-  =  -18dB  (4.9) 


Peaks  at  2.3  Hz  in  the  x  direction,  1.67  Hz,  2.3  Hz,  4  Hz,  4.5  Hz,  and  6  Hz  in  the  y  direction,  and  4 
Hz,  4.5  Hz,  and  6  Hz  in  the  z  direction  have  the  same  spectral  level  in  r(n )  and  record  1207.  Spectral  lev¬ 
els  below  0.5  Hz  in  the  y  and  z  directions  are  also  the  same  in  r(n )  and  record  1207.  Signals  at  these  fre¬ 
quencies  were  not  attenuated  by  the  averaging,  are  therefore  coherent  from  record  to  record,  and  are  like¬ 
wise  attributed  to  the  float’s  resonant  oscillation. 

Also  in  Figure  4.7e,  8  Hz  peaks  in  the  x  and  y  directions  are  about  18  dB  lower  in  r(n)  than  in 
record  1207.  The  18  dB  attenuation  indicates  that  the  8  Hz  signals,  like  the  background  noise,  are 
incoherent  from  record  to  record.  This  further  indicates  that  the  8  Hz  peaks  are  primarily  due  to  ambient 
ocean  noise,  at  least  in  the  x  and  y  directions.  Peaks  at  about  9  Hz  are  also  attenuated  by  about  20  dB  in 
the  x  and  y  directions  and  12  dB  in  the  z  direction,  indicating  that  they  too  originate  in  the  ambient  ocean 
noise  field. 


Table  4.8  -  Spectral  Levels,  Float  3,  Average  of  Records  1201*1260 


Axis 

Frequency  x  y  z 


0.97 

93 

% 

1.5 

92 

- 

- 

1.67 

94 

104 

85 

2.0 

91 

84 

88 

2.3 

100 

96 

86 

4 

79 

72 

76 

4.5 

64 

99 

81 

6 

73 

68 

68 

8 

75 

76 

66 

10 

66 

60 

65 

12 

66 

62 

64 

14 

66 

61 

63 

16 

64 

61 

61 

18 

64 

61 

58 

20 

63 

60 

62 

22 

61 

63 

58 

24 

61 

62 

58 

Of  further  note  in  Figure  4.7e  is  that  many  of  the  2  Hz  harmonics  peaks  are  between  4  and  10  dB 
lower  in  r(n )  than  in  record  1207.  This  indicates  that  their  frequency  or  initial  phase  varies  a  small  amount 
from  record  to  record.  One  plausible  explanation  would  be  that  the  2  Hz  fundamental  frequency  changes 
very  slightly  from  record  to  record.  This  frequency  jitter  could  cause  the  2  Hz  signal  to  remain  essentially 
coherent,  but  the  harmonics  might  change  frequency  over  a  larger  range  proportional  to  the  ratio  of  their 
frequency  to  that  of  the  fundamental.  A  0.01  Hz  change  in  the  2  Hz  signal,  for  example,  would  result  in  a 
0.1  Hz  change  in  the  20  Hz  signal.  Note  that  in  Figure  4.7e  the  2  Hz  signal  is  largely  coherent  from  record 
to  record,  the  12  Hz  is  less  so  (attenuated  about  6  dB  in  the  x  and  y  axes,  no  attenuation  in  the  z  axis),  and 
the  20  Hz  signal  is  largely  incoherent  from  record  to  record  (attenuated  about  9  dB  in  the  x  axis,  18  dB  in 
the  y  axis,  and  8  dB  in  the  z  axis). 

Figures  4.7f  and  4.7g  show  record  1207  after  processing  using  equation  (3.24),  and  Figure  4.7h  com¬ 
pares  the  spectra  of  record  1207  before  and  after  processing.  Processing  reduces  or  removes  spectral  peaks 
attributed  to  resonant  oscillation  at  2.3  Hz,  4  Hz,  and  6  Hz  in  the  x  direction,  1.67  Hz,  2.3  Hz,  4  Hz,  4.5  Hz, 
and  6  Hz  in  the  y  direction,  and  4  Hz,  4.5  Hz,  and  6  Hz  in  the  z  direction.  Processing  also  lowers  spectral 
levels  below  1  Hz  in  the  x  and  y  directions. 

Processing  does  not  affect  the  level  of  spectral  peaks  which  are  attributed  to  the  ambient  noise  field, 
such  as  those  at  8  and  9  Hz.  However,  several  peaks  at  harmonics  of  2  Hz,  which  are  attributed  to  the 
resonant  oscillation,  are  also  unaffected  by  processing,  possibly  because  of  the  aforementioned  fundamen¬ 
tal  frequency  jitter. 

Figures  4.8a  through  4.8h  contain  similar  results  for  record  1607  from  float  3.  The  spectral  com¬ 
parison  plots,  4.8e  and  4.8h,  show  that  processing  reduces  the  level  of  spectral  peaks  attributed  to  resonant 
oscillation  at  2.3  Hz,  4  Hz.  6  Hz,  8  Hz,  10  Hz,  12  Hz,  14  Hz,  16  Hz,  and  22  Hz  in  the  x  direction;  1.67  Hz, 
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2.3  Hz,  4  Hz,  and  4.5  Hz  in  the  y  direction;  and  2.3  Hz,  4  Hz,  4.5  Hz,  6  Hz,  10  Hz,  12  Hz,  16  Hz,  18,  20 
Hz,  and  22  Hz  in  the  z  direction.  Processing  also  reduces  spectral  levels  below  1  Hz  in  the  x  and  y  direc¬ 
tions.  Not  all  of  the  2  Hz  harmonics  attributed  to  resonant  oscillation  are  reduced,  as  in  record  1207.  Sig¬ 
nals  present  in  the  ambient  noise  field  are  not  affected  by  the  processing. 

Figures  4.9  and  4.10  contain  similar  plots  for  records  1207  and  1607  from  float  9,  and  Figures  4.1 1 
and  4.12  contain  results  of  processing  the  same  records  in  float  10.  In  all  cases,  processing  individual 
records  using  a  40  record  average  and  an  estimate  for  a;  removes  or  substantially  reduces  the  major  spec¬ 
tral  lines  attributed  to  resonant  oscillation. 

In  records  1207  and  1607  of  float  10,  processing  attenuates  resonant  oscillation  signals  at  1.5  Hz  in 
the  y  axis  and  4.1  Hz  in  all  three  axes  by  about  16  dB.  However,  the  signals  have  more  than  20  dB  signal 
to  noise  ratio  (SNR)  before  processing  and  are  still  visible  in  the  processed  records. 
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Table  4.1  -  Float  3, 1986  deployment 

y2  ‘2 

"2  M  _.2  *2  •»* 


M 

records 

axis 

i2 

M 

Yr 

It 

Y2 

M 

20 

1201-1220 

X 

.0427 

.002135 

.0239 

.021765 

10.19 

y 

.0613 

.003065 

.0611 

.058035 

18.93 

Z 

.00796 

.000398 

.00453 

.004132 

10.38 

20 

1221-1240 

X 

.0382 

.001910 

.0237 

.021790 

11.40 

y 

.0583 

.002915 

.0599 

.056985 

19.54 

z 

.00588 

.000294 

.00475 

.004456 

15.15 

20 

1601-1620 

X 

.0451 

.002255 

.0281 

.025845 

11.46 

y 

.0863 

.004315 

.0651 

.060785 

14.08 

Z 

.00731 

.000365 

.00336 

.002995 

8.20 

20 

1621-1640 

X 

.0407 

.002035 

.0250 

.022965 

11.28 

y 

.0835 

.004175 

.0605 

.056325 

13.49 

z 

.00703 

.000351 

.00335 

.002999 

8.54 

40 

1201-1240 

X 

.0402 

.001005 

.0233 

.022295 

22.18 

y 

.0580 

.001450 

.0602 

.058750 

40.51 

z 

.00528 

.000132 

.00445 

.004318 

32.71 

40 

1601-1640 

X 

.0429 

.001072 

.0260 

.024928 

23.25 

y 

.0847 

.002117 

.0622 

.060083 

28.38 

z 

.00715 

.000178 

.00327 

.003092 

17.37 

60 

1201-1260 

X 

.0409 

.000681 

.0234 

.022719 

33.36 

y 

.0580 

.000966 

.0604 

.059434 

61.52 

Z 

.00516 

.000086 

.00444 

.004354 

50.62 

60 

1601-1660 

X 

.0437 

.000728 

.0266 

.025872 

35.53 

y 

.0852 

.001420 

.0631 

.061680 

43.43 

z 

.00720 

.000120 

.00328 

.003160 

26.33 
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Table  4.2  -  Float  5, 1986  deployment 


M 

records 

axis 

* 

*2 

r, 

M 

y! 

*2 

Yr 

Y,2 

Yr2 

M 

40 

1001-1040 

X 

.024137 

.000603 

.005064 

.004461 

7.39 

y 

.026275 

.000656 

.0037718 

.0031158 

4.74 

Z 

.0040842 

.000102 

.00035344 

.00025144 

2.46 

40 

1201-1240 

X 

.026102 

.000652 

.0036248 

.0029728 

4.55 

y 

.023268 

.000581 

.0022631 

.0016821 

2.89 

Z 

.0051705 

.000129 

.0003433 

.0002143 

1.66 

40 

1501-1540 

X 

.194700 

.004867 

.010289 

.005422 

1.11 

y 

.072920 

.001823 

.0048317 

.0030087 

1.65 

Z 

.059752 

.001493 

.0017479 

.0002549 

.17 

40 

1701-1740 

X 

.036070 

.000901 

.0031857 

.0022847 

2.53 

y 

.028199 

.000704 

.0022439 

.0015399 

2.18 

z 

.0065993 

.000164 

.00033059 

.00016659 

1.01 

80 

1001-1080 

X 

.025068 

.000313 

.0044149 

.0041019 

13.10 

y 

.028085 

.000351 

.0032843 

.0029333 

8.35 

z 

.0042591 

.000053 

.00030403 

.00025103 

4.73 

80 

1201-1280 

X 

.035478 

.000443 

.0036023 

.0031593 

7.13 

y 

.030123 

.000376 

.0020839 

.0017079 

4.54 

z 

.010001 

.000125 

.00032428 

.00019928 

1.59 

80 

1501-1580 

X 

.181310 

.002266 

.0062007 

.0039347 

1.73 

y 

.085152 

.001064 

.0040127 

.0029487 

2.77 

z 

.052335 

.000654 

.03044 

.029786 

45.54 

80 

1701-1780 

X 

.036614 

.000457 

.0027925 

.0023355 

5.11 

y 

.026466 

.000330 

.0019895 

.0016595 

5.02 

z 

.0060379 

.000075 

.00026261 

.00018761 

2.50 

200 

1001-1200 

X 

.030885 

.000154 

.0040074 

.0038534 

25.02 

y 

.032732 

.000163 

.0026009 

.0024379 

14.95 

Z 

.001521 

.000007 

.00028529 

.00027829 

39.75 

200 

1201-1400 

X 

.043546 

.000217 

.0039816 

.0037646 

17.34 

y 

.086227 

.000431 

.0022422 

.0018112 

4.20 

Z 

.021514 

.000107 

.00035352 

.00024652 

2.30 

200 

1501-1700 

X 

.089202 

.000446 

.003495 

.003049 

6.83 

y 

.087931 

.000439 

.0025187 

.0020797 

4.73 

Z 

.016671 

.000083 

.00039813 

.00031513 

3.79 

Table  4.2  -  Float  5, 1986  deployment 


M 

records 

axis 

Y,2 

1 

M 

y! 

i2 

Y2 

Y2 

M 

200 

1701-1900 

X 

.052966 

.000264 

.0028392 

.0025752 

9.75 

y 

.072233 

.000361 

.0020118 

.0016508 

4.57 

Z 

.021281 

.000106 

.00026349 

.00015749 

1.48 

400 

1001-1400 

X 

.037920 

.000094 

.005064 

.004970 

52.87 

y 

.063715 

.000159 

.0037718 

.0036128 

22.72 

Z 

.012528 

.000031 

.00035344 

.00032244 

10.40 

400 

1501-1900 

X 

.067506 

.000168 

.0029586 

.0027906 

16.61 

y 

.083492 

.000208 

.0020024 

.0017944 

8.62 

z 

.015138 

.000037 

.00027517 

.00023817 

6.43 
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Table  4.. 

M  records  axis 

1201-1220  x 

y 

z 

1221-1240  x 

y 

z 

1601-1620  x 

y 

z 

1621-1640  x 

y 

z 

1201-1240  x 

y 

z 

1601-1640  x 


-  Float  9, 1986  deployment 


0361 

.001805 

.0259 

.024095 

0515 

.002575 

.0439 

.041325 

0107 

.000535 

.00847 

.007935 

0311 

.001555 

.0260 

.024445 

.0313 

.001565 

.0466 

.045035 

,00751 

.000375 

.008265 

0325 

.001625 

.0276 

.025975 

.005165 

.0966 

.091435 

0147 

.000735 

.009465 

.002035 

.0279 

.025865 

iK 

.006000 

.0951 

.089100 

0128 

.000640 

.009360 

0333 

.000832 

.0255 

.024668 

0311 

.000777 

.0444 

.043623 

00825 

.000206 

.00836 

.008154 

0376 

.0273 

.026360 

1151 

.002877 

.0948 

.091923 

0139 

.000347 

.009653 

Table  4.4  -  Float  10, 1986  deployment 


M 

records 

axis 

y? 

Y,2 

M 

a 

*2 

y r 

y* 

M 

20 

1201-1220 

X 

.0624 

.003120 

.0562 

.053080 

17.01 

y 

.190 

.009500 

.177 

.167500 

17.63 

Z 

.0279 

.001395 

.0250 

.023605 

16.92 

20 

1221-1240 

X 

.0479 

.002395 

.0584 

.056005 

23.38 

y 

.1773 

.008865 

.1801 

.171235 

19.31 

Z 

-.0421 

-.002105 

.0262 

.028305 

-13.44 

20 

1601-1620 

X 

.056 

.002800 

.0565 

.053700 

19.17 

y 

.2375 

.011875 

.2000 

.188125 

15.84 

Z 

.04089 

.002044 

.02924 

.027196 

13.30 

20 

1621-1640 

X 

.0675 

.003375 

.0572 

.053825 

15.94 

y 

.2296 

.011480 

.1981 

.186620 

16.25 

z 

.03597 

.001798 

.02937 

.027572 

15.33 

40 

1201-1240 

X 

.05186 

.001296 

.0568 

.055504 

42.82 

y 

.1823 

.004557 

.1775 

.172943 

37.95 

z 

.00093 

.000023 

.0251 

.025077 

1090.30 

40 

1601-1640 

X 

.0611 

.001527 

.0564 

.054873 

35.93 

y 

.2322 

.005805 

.1985 

.192695 

33.19 

Z 

.03734 

.000933 

.02917 

.028237 

30.26 
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Figure  1.2 
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